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Abstract 

After being destroyed by a binary supermassive black hole, a stellar density cusp can regrow at the center 
of a galaxy via energy exchange between stars moving in the gravitational field of the single, coalesced hole. 
We illustrate this process via high-accuracy A^-body simulations. Regeneration requires roughly one relaxation 
time and the new cusp extends to a distance of roughly one-fifth the black hole's influence radius, with density 
p ~ r~ 7 / 4 ; the mass in the cusp is of order 10% the mass of the black hole. Growth of the cusp is preceded 
by a stage in which the stellar velocity dispersion evolves toward isotropy and away from the tangentially- 
anisotropic state induced by the binary. We show that density profiles similar to those observed at the center of 
the Milky Way and M32 can regenerate themselves in several Gyr following infall of a second black hole; the 
presence of density cusps at the centers of these galaxies can therefore not be used to infer that no merger has 
occurred. We argue that p ~ r~ 7 / 4 density cusps are ubiquitous in stellar spheroids fainter than My w —18.5 
that contain supermassive black holes, but the cusps have not been detected outside of the Local Group since 
their angular sizes are less than ~ 0.1". We show that the presence of a cusp implies a lower limit of ~ 10~ 4 
yr -1 on the rate of stellar tidal disruptions, and discuss the consequences of the cusps for gravitational lensing 
and the distribution of dark matter on sub-parsec scales. 
Subject headings: 

1. INTRODUCTION 

Mass distributions near the centers of early-type galaxies 
are well described as power laws, p ~ r~V, with indices y 
that change gradually with radius. At their innermost resolved 
radii, most galaxies have 0.5 ^ y ^ 2., with the steeper slopes 
characteristic of fainter galaxies (tFerrarese et al.ll2005l) . If a 
supermassive black hole is present, the orbits of stars will be 
strongly influenced at distances less than ~ rj, = GM./a 2 ss 
10 pc (M./10 8 M@)(a/200km s" 1 )" 2 , the black hole's grav- 
itational influence radius. Most galaxies are spatially unre- 
solved on these small scales; two clear exceptions are the 
nucleus of the Milky Way, for which nu mber counts extend 
inward to ~ 0.002^ JGenzel et alJl2003l). and M32, w hich is 
resolved down to a radius of ~ 0.2r/, ( Lauer et al. 1998). Both 
galaxies exhibit steep density slopes, 1.5, at r r/,. Out- 
side of the Local Group, only giant ellipticals have sufficiently 
large black holes that r/, can be resolved; the nuclear luminos- 
ity profiles in these galaxies are also power laws but very flat, 
Y^l. 

Many distributions of stars are possible around a black hole, 
but under two circumstances, the stellar distribution at r & r n 
is predictable. (1) If the black hole has been present for a 
time longer than T r , the relaxation time in the nucleus, ex- 
change of energy between stars will drive the stellar distri- 
bution toward a collisional steady state; assuming a single 
stellar mass and ignoring physical collisions between stars, 
this s teady state has p ~ r~ 7 / 4 al r <, r/, (Ba hcall & Wolfl 
119761) . (2) If the nucleus formed via the merger of two galax- 
ies each with its own supermassive black hole, the two black 
holes will displace of order their combi ned mass in the pro- 
cess o f forming a tightly -bound pair ( Milos avTievic & Merittl 
1200 ll) . producing a low-density core. The first mechanism 
may be responsible for the steep density profiles observed at 
the centers of the Milky Way and M32, since both galax- 
ies have central relaxation times of order 10 9 yr and both 
are near enough that linear scales of order r/, are well re- 



solved. The second mechanism may explain the very flat cen- 
tral profiles of luminous E galaxi es ( Milosavlievic et alJ2002t 
iRavindranath. Ho"&^ Filippenko 2002); the central relaxation 
times of these galaxies are much longer than 10 10 yr and the 
stellar distribution would be expected to remain nearly un- 
changed after the two black holes had coalesced into one. 

In this paper we point out that both outcomes are possi- 
ble. A galaxy may form via mergers, but at the same time, its 
central relaxation time following the merger may be shorter 
than 10 10 yr. In this circumstance, the cusp of stars around 
the black hole is first destroyed by the massive binary, then is 
regenerated via encounters between stars in the gravitational 
field of the single, coalesced hole. The result is a steep inner 
density profile in a galaxy that had previously experienced 
the scouring effects of a massive binary. To the extent that all 
stellar spheroids experienced mergers - if only in the distant 
past - this picture is probably generic, applying even to small 
dense systems like M32 and to the bulges of spiral galax- 
ies like the Milky Way. Understanding the conditions under 
which a previously-destroyed density cusp can spontaneously 
regenerate is crucial if one wishes to interpret the present-day 
luminosity profiles of galaxies as fossil relics of their merger 
histories llVolonteri et all2003l) . 

We use A^-body simulations (§2) to follow first the destruc- 
tion (§3) and then the spontaneous regeneration (§4) of den- 
sity cusps around black holes. The two most important free 
parameters in this problem are the mass ratio q = MijM\ of 
the binary black hole, and the slope y of the initial density 
cusp surrounding the larger hole. We present results for sev- 
eral combinations of q and y (Table 1). Our conclusion is that 
collisional, Bahcall-Wolf density cusps should be ubiquitous 
in stellar spheroids fainter than My w — 1 8 .5 that contain mas- 
sive black holes, essentially regardless of their merger histo- 
ries. However these cusps have gone undetected in galaxies 
outside the Local Group because they are unresolved. In §5 
we discuss a number of consequences of the presence of the 
cusps. 
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TABLE 1 

Parameters of the A'-body integrations 



Run 


1 


M 2 /Mi 


'in 


'ftp 


TJn ) 


ah 


ri 


TJri) 


T r (0 2r[ ) 


T 

sop 


r'! 
h 


1 


0.5 


0.5 


0.264 


0.326 


1170. 


0.0181 


0.39 


1420. 


620 


79. 


0.38 


2 


0.5 


0.25 


0.264 


0.296 


1010. 


0.0119 


0.35 


1210. 


420 


48. 


0.34 


3 


0.5 


0.1 


0.264 


0.278 


916. 


0.00573 


0.32 


1070. 


390 


22. 


0.31 


4 


1.0 


0.5 


0.165 


0.210 


599. 


0.0116 


0.28 


870. 


300 


48. 


0.27 


5 


1.0 


0.25 


0.165 


0.188 


499. 


0.00751 


0.23 


640. 


270 


26. 


0.23 


6 


1.0 


0.1 


0.165 


0.174 


441. 


0.00360 


0.20 


520. 


220 


11. 


0.20 


7 


1.5 


0.5 


0.0795 


0.107 


217. 


0.00594 


0.17 


420. 


160 


23. 


0.17 


8 


1.5 


0.25 


0.0795 


0.093 


170. 


0.00364 


0.13 


260. 


110 


10. 


0.13 


9 


1.5 


0.1 


0.0795 


0.085 


144. 


0.00176 


0.11 


200. 


75 


4.1 


0.11 



2. MODELS AND METHODS 

We started by constructing Monte-Carlo realizations of 
steady-state galaxy models having Dehnen's (1993) density 
law, with an additional, central point mass representing a 
black hole. The Dehnen-model density follows p(r) °< r~T 
at small radii, and the isotropic phase-space distribution func- 
tion that reproduces Dehnen's p(r) in the presence of a cen- 
tral point mass is non-negative for all y > 0.5; hence y = 0.5 
is the flattest central profile that can be adopted if the initial 
conditions are to represent an isotropic, steady state. We con- 
sidered initial models withy= (0.5,1.0,1.5). The massMi of 
the central "black hole" was always 0.01, in units where the 
total mass in stars M ga i was one; the Dehnen scale length ro 
and the gravitational constant G were also unity. The A^-body 
models so constructed were in a precise steady state at time 
zero. 

Destruction of the cusp was achieved by introducing a sec- 
ond "black hole" into this model, which spiralled into the cen- 
ter, forming a binary with the first (more massive) hole and 
displacing stars. Three values were used for the mass of the 
smaller hole: M 2 /M x = q = (0.5,0.25,0.1). The smaller hole 
was placed initially at a distance 1.6 from the center, with a 
velocity roughly 1/2 times the circular velocity at that radius; 
a non-circular orbit was chosen in order to speed up the orbital 
decay. 

After the orbit of the smaller black hole had decayed via dy- 
namical friction against the stars, it formed a tight binary with 
the more massive hole, with a relative orbit close to circular. 
An estimate of the semi-major axis a/, at which the binary first 
becomes "hard" is a/, = G/j/4o 2 where p = M\M\j(M\ +M2) 
is the reduced mass. The precise meaning of "hard" is debat- 
able; the definition just given defines a "hard" binary as one 
whose binding energy per unit mass, \E\/(M\ +M2), exceeds 
2a 2 . While simple, this definition contains the ill-defined 
quantity o, which is a steep function of positi on near the black 
hole(s). We followed lMerritt & Wangl l l2005i) and used the al- 
ternative definition 

_ n r hn _ q r hn 

h ~Mi+M 2 4 (l+qf 4 ' K ) 

with o, 12 the gravitational influence radius defined below. In 
practice, a/, so defined was found to be roughly (within a fac- 
tor <~ 2) the value of the semi-major axis at which the binary 
hardening rate (d/dt)(\ /a) first became approximately con- 
stant. 

Decay was allowed to continue until the binary semi-major 
axis had reached a value of a/,/5 « qr^j '20. At this point, 



the two black holes were replaced by a single particle of mass 
M\2=M\ +M2, with position and velocity given by the center 
of mass of the binary. The A^-body integration was then con- 
tinued for a time roughly equal to the relaxation time r r (r^) 
defined below. The most suitable time at which to merge the 
two black holes was not clear a priori; our choice is of order 
the separation at which gravitation al-wave emission would 
induc e coalescence in ~ 10 10 yr jMerritt & Milosa vTievid 
2005), but in fact we expect that other processes like inter- 
action of the binary with ambient gas may drive the final coa- 
lescence in real galaxies, and it is not clear at what separation 
these processes are likely to dominate the evolution. 

The ratio (M\ + Mi) / 'M ga i ~ 0.01 in our models is roughly 
a factor ten larger than the rati o of black hole mass to 
galaxy mass in real spheroids jMerritt & Ferraresel 1200 It 
Marc onT& Huntll2003l) . This is acceptable as long as we are 
careful to present masses and radii in units scaled to M\ +M2 
when making comparisons with real galaxies. 

Table 1 gives a number of parameters associated with the 
Af-body integrations. The gravitational influence radius r/, was 
defined as the radius containing a mass in stars equal to twice 
the mass M. of the central black hole. This definition, while 
superior to GM./o 2 , is somewhat ambiguous in our A^-body 
models, given that the effective mass of the central object, 
and the distribution of the stars, both change with time. We 
accordingly defined four different influence radii. (1) At the 
start of the integrations, the larger black hole, of mass M\ , 
was located at the center of the galaxy. Its influence radius r/ H 
was computed by setting M. = M\ and using the t = stellar 
distribution. (2) After the smaller black hole has fallen in to 
a distance <, r^ from the larger hole, the appropriate value 
of M. becomes M\ +M2. We defined the associated influence 
radius to be r/, 12 , which we computed ignoring the changes 
that had occurred in the stellar distribution since t = 0. (3) The 
third influence radius, ri, was computed by setting M. = Mi + 
M2, but this mass was compared with the stellar distribution at 
the end of the binary evolution phase, after the phase of cusp 
destruction. (4) Finally, r" h is the influence radius at the end 
of the second phase of integration, after a Bahcall-Wolf cusp 
has formed around the single black hole. As Table 1 shows, 
r^' is only slightly smaller than r' h since the regenerated cusp 
contains a mass that is small compared with M.. 

The relaxation times T r in Table 1 were computed from 
the standard expression (eq. 2-62 of Spitzer 1987), setting 
InA = ln(o 2 r/,i2/2Gm*), with a the ID stellar velocity dis- 
persion at r/,^ and m 4 = A^ 1 the mass of an A^-body particle. 
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FIG. 1. — Cusp destruction. Solid lines are stellar density profiles just before the two "black holes" were combined into one. Dashed lines show the initial 
models, q = M2/M1 is the binary mass ratio. 



This definition of A is equivalent to equating b max , the max- 
imum impact pa rameter for encounters in Chandrasekhar's 
theory, with r/, llPreto. Merritt & Spurzeml 120041) . Table 1 
gives values of T r evaluated at two of the four influence radii 
defined above. T(rh n ) was computed using the structural pa- 
rameters of the initial galaxy model, while T)-(r^) was com- 
puted using estimates like those in Figures [2 and [3] of the 
stellar density and velocity dispersion in the evolved models. 
We also give T r (0.2rh'); the motivation for this is given be- 
low. The final time scale in Table 1, T gap , is defined below 
(equation and is an estimate of the time required for the 
angular-momentum gap created by the binary to be refilled; 
Tgap varied between ~ T r (ri, )/20 and ~ T r (r^)/5Q. 

The power-law cusps in our initial models were motivated 
by the approximately power-law dependence of luminosity 
density on radius o bserved near th e centers of many early- 
type galaxies (Ferra rese et alJl2005l) . Since the observations 
often do not resolve r^, the stellar distribution at r < r/, in 
some galaxies might be different than the inward extrapola- 
tion of the power laws that are fit to larger radii. For instance, 
galaxies with sufficiently short relaxation times are expected 
to have p ~ r~ 7 / 4 density cusps like the ones that form in 
our jV-body models at late times. Other galaxies may have 
compact stellar nuclei (ICote et al.l l2005). We did not include 
such dense features in our initial models: first, because the 
associated mass would have been small compared with the 
mass removed by the binary; and second, because doing so 
would have more than doubled the computational effort due 
to the short time steps required for stars initially near the black 
holes. 

All N-body integrations used N = 0.12 x 10 6 par- 
ticles and were carried out on a GRAPE-6 special- 
pu rpose computer. The jV-body integrator is described 
in IMerritt. Mikkola & SzeiH (120061). This algorithm is an 
adaptation of NB0DY1 lAarsethl d!999h to the GRAPE-6; it 
uses a fourth-order Hermite integratio n scheme wi th indi- 
vidual, adaptive, block time steps (lAarsethl 12003). For 
the majority of the particles, the forces and force deriva- 
tives were calculated via a direct-summation scheme on the 
GRAPE-6, using the particle advancement scheme described 



in iBerczik. Merritt & Spurzeml (12005) with an accuracy pa- 
rameter of T| = 0.01 and zero softening. Close encounters be- 
tween the black holes, and between black holes and stars, re- 
quire prohibitively small time steps in such a scheme and were 
regularized using th e chain regularization rou tine of Mikkola 
and Aarseth (Mikkola & Aars efhll990lll993l) . A detailed de- 
scription of the chain algor ithm, inc luding t he results of per - 
formance tests, are given in Merritt, Mikkola & Szell (2006). 

Our initial conditions (one black hole at the center, a smaller 
black hole orbiting about it) are not as realistic as in simula- 
tions that follow both merging galaxies from the start (e.g. 
Milos avlievic & Merittl J2001), M erritt et all ll2002l) 1. but are 
superior to simulations that drop one or two black holes 
into a pre-existing galaxy th at contains no black hole (e. g. 
lOuinlan & Hernauisd dl997l) : Nak ano & Makinol ( Il999albl) '). 
We ignore the radiation recoil that would accompany the final 
coalescence of the two black holes, displacing the remnant 
hole temporarily from its central location and increasing the 
size o f the core ( Merritt et alJl2004bt iBovlan-Kolchin et alJ 
2004). We also ignore processes like loss of stars into the 
black hole(s), stellar tidal disruptions, and stellar collsions, 
all of which might affect the form of the final density profile. 

3. CUSP DESTRUCTION 

After formation of a hard binary at a a/,, the bi- 
nary's binding energy continues to increase as the two 
massi ve particles eject stars via the gravitational sling- 
shot (ISaslaw. Valtonen & Aarsethlll974tlMikkola & Valtonenl 
1992; Ouinlan 1996). As in a number of recent studies 
(IBerczik. Merritt & Spurzemll2005l IMerritt. Mikkola & Szell 
2006), we found that the hardening rates were nearly in- 
dependent of time, s = (d / dt)(\ / a) « const., for a < ah. 
The values of s are not of particular interest here, since N- 
body harden ing rates are known to depend stro ngly on par- 
ticle number! Berczi kTMerritt & Spurzemll20 05). and we do 
not expect these simulations to be in the "empty loss cone" 
regim e that characterizes bin ary evolution in real galaxies 
(Milosavlievic & Meritt 2003). However the "damage" done 
by the binary to the pre-existing stellar distribution is expected 
on energetic grounds to be a function of essen tially in- 
dependent of the rate of hardening and hence of N (|Quinlan 
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FIG. 3. — Stellar velocity dispersions at the end of the binary phase, in the integrations with M%jM\ = 0.5. Solid lines: o,-; dashed lines: o~ f . The binary creates 
a tangentially-biased velocity distribution near the center by preferentially ejecting stars on radial orbits. Arrows indicate the black hole influence radius r' h . 



1996), and this expe ctation has recently been confi rmed in 
iV-body experiments (Merritt, Mik kola & Szell20 06). As de- 
scribed above, the binary phase of the integrations was termi- 
nated when a = ah/5. 

Figure^shows density profiles of the A^-body models at the 
end of the binary phase. Radii of the "star" particles were 
computed relative to the center of mass of the binary, and 
estimates of p(r) were constructed using the nonparame trie 
kernel estimator of Merritt & Tremblav ( 1994). The damage 
done by the inf ailing black hole can be seen to increase with 
its mass for y = 0.5 and 1, and extends out to a radius ~ r/,. 
The y = 0.5 cusp is converted into an approximately constant- 
density core, while the y = 1 and y = 1 .5 cusps are "softened," 
to power laws of index ~ 0.5 (y — 1) and ~ 1.0 (y — 1.5). 

A standard measure of the damage done by an infalling 
black hole to a pre-exising density cusp is the "mass deficit" 
Mdef> defined as the decrease in the ce ntral mass w ithin a 
spher e that contains the affected region ( Milosavlievic et al. 
2002). The mass deficit is potentially observable 



( Milosavlie vic et alJ 120021 iRavindran ath. Ho & Filinpenko 
2002; Graham 2004), assuming that one can guess the pre- 
existing density profile, and its value is an index of the cumu - 
lative effect of mergers on the galaxy ( Volonteri et al. 2003). 
Since there are few results in the literature on the sizes of 
mass deficits generated by large-mass-ratio inspirals, we show 
in Figure |3 mass deficits for these integrations, expressed in 
terms of a^/ait); we extended some of the integrations be- 
yond a = a/,/5 in order to further elucidate this dependence. 
Figure|2] shows that - for a given degree of binary hardness - 
the mass deficit is much better predicted by M\ or M\i than 
by M%. In other words, the damage done by the smaller black 
hole is rougly proportional to the mass of the binary, at a given 
value of af,/ a. 

This is a reasonable result. Equating the increase in the 
binary's binding energy with the energy lost to an ejected star 
gives 

SflJLi^ (!)_(,, (2) 
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FIG . 4 . — Mean velocity anisotropy of the stars in a sphere of radius ri around the single black hole, as a function of time after the massive binary was combined 
into a single particle. Shown are the integrations with MijMy = 0.5. Arrows indicate the time T gap defined in the text (equation^ Table 1). 
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FIG. 5. — Cusp regeneration. Plots show the stellar density profile around the single black hole, at post-coalescence times of approximately (0,0.1,0.5, l)r,-(r;,) 
(y = 0.5), (0,0.2, l)T r (r/,) (y = 1.0), and (0,0.2,0.5)r r (r;,) (y = 1.5). Dotted lines show the density profiles at the start of the binary integrations, when only one 
black hole was located at the center. Dashed lines have logarithmic slope of —1.75. Arrows indicate ri, the influence radius of the single black hole just after 
coalescence of the binary (Table 1). 
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FIG. 6. — Final density profiles after growth of the cusp. Black lines had 
Y = 0.5 for the pre-binary initial conditions, red lines had y = 1.0, and blue 
line had y= 1.5. Dashed lines show the initial models. 



with Ei and Ef the initial and final kinetic energy per unit mass 
of the ejected star and ?>M e j its mass; jj is the reduced mass of 
the binary. This can be written 



dM, 



<-■) 



l-t 



V, 



bin 



<iln(l/a) 2E f -Ei 



(3) 



where Vun — y/G(M\ +M2)/a is the relative velocity be- 
tween the components of the binary (assuming a circular or- 
bit). Hills (1983) finds from scattering experiments that, for a 
hard binary in the limit M2 <C M\, Ef — Ei for a zero-impact- 
parameter encounter is ~ lAG/j/a. This implies 

dM e i 

^h k036{Mi+M2) (4) 

and the total mass ejected by the binary in hardening from 

a m a/, to a/ is 



AM 



Mi +M 2 



i0 .361n [ a ± 

\ a f 



(5) 



i.e., the mass ejected in reaching a given degree of hardening 
an /a is proportional to Mi +M2. Quinlan (1996) used scat- 
tering experiments to compute J in the relation 



dM e 



d\n{l/d 



= /(Mi+M 2 ) 



(6) 



as a function of M2/M1, assuming a Maxwellian distribution 
of velocities at infinity. He defined "ejected" stars to be those 
with final velocities exceeding max{ 1 .5Vo, V3o}; here "final" 
means after the star has moved far from the binary (ignoring 
any velocity changes that result from the galactic potential), 
Vb is the velocity before the encounter, and a is the velocity 
dispersion far from the massive binary. Quinlan's definition 
is a reasonable one if M e j is to be equated with Md e f, and 
Quinlan found that J is nearly independent of M%/M\ for a 
hard binary, decreasing from J « 1 for M2 = Mi to J w 0.5 for 



M? = Mi/256 (his Figure 5). These /-values are consistent 
with Figure 12 

Expressions like (0 seem counter-intuitive, since they seem 
to imply an ejected mass that is large even for small M2. How- 
ever for small M2, a/, °< Mi/M\ (equation[0, so for a given a/, 
the ejected mass is predicted to scale as ~ lnM2. The depen- 
dence of Md e f on binary mass ratio is a topic that deserves fur- 
ther study since the central structure of bright elliptical galax- 
ies is ex pected to reflect the cu mulative effect of many minor 
mergers ( Volont eri et al.l2 003). 

We do observe an additional dependence of M c / e f on both 7 
and q at given ah/ a (Figure|2). However the definition of a/, is 
somewhat arbitrary (§2) and we do not pursue that additional 
dependence here. 

The binary preferentially ejects stars on eccentric orbits. 
The result is an induced anisotropy in the stellar velocity dis- 
tribution at r <, r/,. Figure [3] illustrates this, for the integra- 
tions with the most massive secondaries (M2/M1 = 0.5). The 
induced anisotropy is mild, similar to what has been observ ed 
in other A-body studies (e.g Milosavljevic feMerittl(l2"00"TV ). 

4. CUSP REGENERATION 

When the separation between the two components of the bi- 
nary had dropped to a = ah/5, the two massive particles were 
combined into one and placed at the center of mass of the bi- 
nary. The integrations were then continued until an additional 
time of roughly r r (rl) had elapsed (Table 1). During this sec- 
ond phase of the integrations, the stellar distribution around 
the single black hole evolved due to star-star encounters. The 
initial effect of encounters was to refill the phase space gap 
created when the binary ejected stars with pericenters r p <, a^ 
The refilling time is approximately 



a h 



T r (r h ) : 



1 



4(1+?) 



-T r (r h ) 



(7) 



(Merr itt & Waiig|l2005l) . Table 1 shows that T gap varies be- 
tween ~ 10(y= 1-0,9 = 0.1) and -80 (7=0.5,^ = 0.5) in 
our units, or between — 0.02r r (r/j) and — 0.05r r (7/ 7 ). Fig- 
ure [4] verifies that the initially anisotropic velocity distribu- 
tion at r < r/ 7 becomes isotropic in a time of order T gap . In 
luminous galaxies, eve n T gap can exceed a galaxy lifetime 
(Merritt & Wans 2005), and such galaxies will experience a 
greatly reduced rate of interaction of stars with the black hole 
compared to galaxies with a full loss cone. 

On the longer time scale of — T r (rh), energy exchange 
between stars modifies their radial distribution, eventually 
reaching the zero-flux Bahcall-Wolf (1976) solution p — 
r -7/4 -p n j s j s snown m Figure |3 for the 7 = 0.4 and 7=1 
models. (Results for 7= 1.5 are shown in Figures [6] and Q 
and are discussed separately below.) The new cusp attains the 
p — r~ 7 / 4 form inside of — 0.2r^. In all of the integrations, 
the logarithmic slope of the final cusp is steeper than that of 
the initial (pre-binary) galaxy at radii r ^ 0.5r^, although the 
final density remains below the initial density at r w r J h . 

Figure[6]shows the final density profiles for all seven mod- 
els. There is remarkably little dependence of the final result 
on the mass ratio of the binary. The reason can be seen by 
comparing Figures^and[6] while the different mass ratios do 
yield different density profiles at the end of the binary phase, 
most of these differences are at r ^ 0.3r/ 7 , and this region is 
efficiently refilled by the new cusp. 

The final density profiles in Figure [6] can be well approxi- 
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mated as broken power laws at radii r <J r/,: 

-7/4 



pM = p(n>) 
=p(n>) 



ro < r < r/,. 



(8) 



The "break" radius ro is ~ 0.15^" in the models with initial 
cusp slope y = 0.5, - 0.20r£ for y = 1.0, and - 0.25^' for 
y = 1.5. In this approximation to p(r), the mass in the cusp, 
M cusp = M{r < r ), is given by 



M, 



2M. 



cusp 



1+0.75 



(rh/ro) -1 



(9) 



0.16 



Thus M cusp /M. re 0.03 (y= 0.5), ~ 0.12 (y= 1.0), and > 
(y = 1.5), confirming thatM £! „ p is small compared with M.. 
One consequence is that rl', the influence radius at the final 
time step, differs only slightly from rl, the influence radius at 
the end of the binary phase (Table 1). 

Finally, we consider how much time is required to regrow 
the cusps. As noted above, all of the post-binary integrations 
were carried out until an elapsed time of ~ Tr(r^), and this 
was also roughly the time required for the cusp to reach its 
steady-state form. However the cusp extends only to a radius 
of ~ 0.2rJ ; , and a more relevant estimate of the relaxation time 
is probably r r (0.2r^). Estimates of this time are given in Ta- 
ble 1 . The time to fully grow the cusp is approximately 2 — 3 
times the relaxation time at 0.2rJ r As Figure|5]shows, consid- 
erable regeneration takes place even in a much shorter time, 
hence we might predict the presence of steep central profiles 
even in galaxies substantially younger than either T r {rh) or 
(2 — 3)r, (0.2r/,). We return to this issue below. 

5. IMPLICATIONS 

5.1. Cusp Regeneration in Local Group Galaxies 

The relaxation time in the Milky Way nucleus is less than 
10 10 yr and a number of authors have suggested that the stellar 
cluster arou nd the supermassive black hole might be collision- 
ally relaxed ( Alexander! 19991: iGenzel et alJ2 003). The stellar 
mass density, based on number counts that extend down to 
~ 0.005 pc, is 



p(r): 



1.2 x 10 6 M o pc~ 3 



(10) 



v 0.38pc, 

JGenzel et alJ2003l) . with a pa 2.0 at r > 0.38 pc and a pa 1 .4 
at r < 0.38 pc. T he black hole mass is 3.7 ±0.2 x 10 6 M e 
JGhez et alJl2005l) . The implied influence radius, defined as 
the radius containing a mass in stars that is twice the black 
hole mass, is r/, pa 88" pa 3.4 pc. 

We used equation JlQi and the measured value of M, to 
compute the stellar velocity dispersion profile in the Milky 
Way nucleus, assuming a spherical velocity ellipsoid. The re- 
laxation time T r was then calculated as in the TV-body models, 
i.e. 



Tr(r) 



0.338a(r) 



p(r)m*G 2 In A 
1.69 x 10 9 yrx 

Q(r) W p(r) 
100 km s- 1 / V 10".V/. 



0.7M, 



(11a) 
(lib) 



with 111+ = O.IMq and A = a 2 (r/,)r/ ! /2Gm*. The relaxation 
time at r = r/, was found to be 5.4 x 10 10 yr, dropping to ~ 
6 x 10 9 yr at 0.2r/, and ~ 3.5 x 10 9 yr at ~ O.lr/,; at smaller 
radii, T r increases very slowly with decreasing r. 

As noted above, cusp regeneration in the TV-body models 
required a time of ~ ir r (r/,), or ~ 2— 3r,(0.2r/j). However 
these two timescales are substantially different in the Milky 
Way, suggesting that the structure of the Milky Way nucleus 
differs in some respect from that of the TV-body models in the 
region 0.2r/, & r&r h . 

Figure makes the comparison. We show the evolution 
of p(r) and T.(R) (the latter is the projected density) in the 
post-binary integration with parameters (y = 1.5, q = 0.5); 
the right-hand frame also shows the measured surface den- 
sit y of stars at t he Ga lactic center, from the number counts 
of IGenzel et alJ d2003h . with arbitrary vertical normalization. 
The TV-body model successfully reproduces the slope of the 
observed stellar cusp at late times, but the TV-body profile is 
less steep than the observed counts at r > 0.2r/j, i.e., outside 
of the cusp. In the Milky Way, an approximately power-law 
dependence of p on r extends well beyond r/,, while the TV- 
body models exhibit a shallower slope at these radii. Better 
correspondence between TV-body model and data could pre- 
sumably have been achieved by modifying the initial model, 
or by repeating the integrations with smaller M. (which would 
have mandated a larger TV in order to resolve the smaller 
cusp); the unphysically large value of M, adopted here (M, = 
O.QISMgui) implies an influence radius that is comparable to 
the Dehnen-model scale length. 

In light of these ambiguities, it is reasonable to fix the phys- 
ical unit of time by comparing the TV-body model to the Milky 
Way at ~ 0.2r/,; this is roughly the outer radius of the cusp, 
and Figure0shows that the model fits the data well out to this 
radius. We find that T, (0.2r/,) pa 160 in TV-body units, while 
in the Milky Way, T r {Q.2r h ) pa 6 x 10 9 yr. The integration 
time of 300 in Figure then corresponds to a physical time 
of ~ 1.1 x 10 10 yr. In fact, the TV-body cusp attains nearly its 
final form by a time of ~ 200, or ~ 8 x 10 yr. We conclude 
that the distribution of stars near the Milky Way black hole is 
consistent with a merger having occurred at a time ^ 8 Gyr 
in the past. 

In fact, the probability is ^68% that a galaxy with a dark- 
matter halo as large as that of the Milky Way has experienced 
a major merger (mass ratio of at l east 4:1) since a re dshift of 
z = 2, i.e., in the last ~ 10 Gyr jMerritt et al]l2002l) . If this 
occurred, our results suggest the interesting possibility that 
the Milky Way cusp may still be evolving toward its steady- 
state form. Evidence for this might be sought in the form of 
a flattening the density profile at very small radii, r <, 0. lr/, pa 
10". 

The stellar density near the center of M32 is 



p(r) Pa 2.2 x lO 5 M p C - 3 [j^j 



a pa 1.5 



(12) 



jLauer et alJ 19981) . very similar to that of the Milky Way. Un- 
fortunately the estimation of M, in M32 suffers from the de- 
generacy inherent in orbital modelling of axisymmetric sys- 
tems and the black hole mass could lie anywhere in the 
range 1 .5 x 10 6 M W <j M. <j 5 x 10 6 M W with equal likelihood 
( Vallu rLMerritt & Emselleml20 04). Applying the A/. - o re- 
lation JFerrarese & Merrittl 120001) gives M. pa 2.0 x 10 6 M e 
which we adopt here. The black hole influence radius be- 
comes ~ 1.7 pc pa 0.5" and T r {r h ) pa 2.2 x 10 10 yr. This is 
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FIG. 7. — Evolution of the spatial (a) and projected (b) density profiles in the simulation that most closely resembles the Milky Way nucleus. The initial model 
(dotted line) had a p ~ r~' 5 density cusp and a central black hole of mass Mi = 0.01, in units where the total stellar mass of the model is 1. Blue (lower) line 
shows the density after infall and capture of a second black hole of mass Mi = 0. \M\ . Green and red lines show the evolving density at times (40, 100) after 
coalesence; based on the discussion in the text, the corresponding physical times are roughly (3,8) x 10 9 yr. Symbols in (b) show the observed surface density of 
stars near the Galactic Center, from Genzel et al. (2003); see Figure 7 from that paper for symbol definitions. The vertical normalization of the symbols was set 
by eye in order to give a good match to the iV-body data at R < 0.2r;, . Dashed lines have logarithmic slopes of — 1 .75 (a) and —0.75 (b). 



somewhat smaller than the estimate of T r (r^) in the Galactic 
nucleus, suggesting that a collisional cusp could easily have 
been regenerated in M32. However, the cusp should only ex- 
tend outward to ~ 0.2r/, ss 0.1" at the distance of M32, barely 
resolvable even with HST. Depending on the exact value of 
M., the collisional cusp in M32 may or may not have been 
resolved. 

Both the mass of the black hole in M31, and the dynamical 
state of its nucleus, are less certain than in the Milky Way 
and M32 due to M31's complex morphology. However most 
estimates of the relaxation time at the center of P2, the nuclear 
component believed to contain the black hole, are of order 
10 11 yr (e.g. Lauer et al. 1998). M33 has a very short central 
relaxation time but no dynamical signature of a black hole 
JVallurietal.120051). 

5.2. Cusp Regeneration in Galaxies Beyond the Local Group 

Which galaxies beyond the Local Group should contain 
Bahcall-Wolf cusps? Figure[8]shows estimate s of T r (rh) in the 
sample of early-type galaxies modelled bv IWang & Merritil 
(2004); these are a subset of the Magorrian et al.llll99 8) galax- 
ies. In this plot, filled symbols denote galaxies observed with 
angular resolution & j, s < B r ,, and the size of the symbol is 
proportional to log(8, 7i /9 6 J ). (In galaxies where r/, is not 
resolved, we note that the estimates of T r are particularly un- 
certain, since they depend on an inward extrapolation of lu- 
minosity profiles measured at r > r/,.) 

Figure [8] reveals that the brightest spheroids, My <, —20, 
have central relaxation times that always greatly exceed 10 10 
yr. In these galaxies, a low-density core created by a bi- 
nary supermassive black hole would persist for the age of 
the universe. However r,-(r/,) drops with decreasing luminos- 



ity, falling below 10 10 yr for My ^ - 18. The Milky Way 
bulge falls on the relation defined by the more distant galaxies, 
which is reassuring given the uncertainties in its luminosity. 
(We adopte d a bulge blue abs olute magnitude of Mb = —17.6 
from lMarconi & Hurl J2Q03) and assumed My = M B — 0.9.) 
However M32 appears to be shifted from the relation defined 
by the other galaxies, as if it is the dense core of a once much 
bright er galaxy. This possibility has often been raise d in the 
past {king 1962t lFaberll973l:lNieto & Prugnielll987l) . 

Figure [8] combined with the arguments in §5.1, suggests 
that spheroids fainter than My ps —18.5 are dynamically old 
enough to contain Bahcall-Wolf cusps. 

However outside of the Local Group, Figure [8] also sug- 
gests that such cusps are unlikely to be resolved. We can 
check this prediction by relating the angular size of the cusp 
to black hole mass via the M. — o relation, M% ps 1 .66O200 6 
JFerrarese & Fordll2005h . with M g = M./1O 8 M and a 2 oo 
O 200 km s -1 , and (temporarily) redefining r/, as GM m /a 2 ps 
11.2M8a2 2 pc. Then 

r h ps 2.76MjP 9 pc. (13) 

Taking for the outer radius of the cusp ro ps 0.2r/, (equation 
[8J1, its angular size becomes 

e «0.57"M 8 °- 59 D M ; x . (14) 

with Dm P c the distance to the galaxy in Mpc. Thus, at the dis- 
tance of the Virgo cluster, 0o > 0. 1" implies M. <; 6 x 1O 8 M ; 
however such massive black holes would almost certainly sit 
in galaxies with central relaxation times longer than 10 10 yr 
(Figure [8jl and a cusp would not have formed. If we assume 
that black holes like the ones in the Milky Way and M32 (i.e. 
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FIG. 8. — Estimates of the relaxation time at the black hole's influence ra- 
dius, o,, in the sample of early-type galaxies modelled by Wans & Merritt 
(2004). Black hole mas ses were computed from the M, — a relation 
i Merritt & Ferrarese 2001), except in the case of the Milky Way, for which 
M. = 3.7 x 10 6 M o was assumed iGhez et al. 2005). The stellar mass was 
set equal to O.7M when computing T r . Horizonal axis is absolute visual 
magnitude of the galaxy or, in the case of the Milky Way, the stellar bulge. 
The size of the symbols is proportional to log| (9, Vi /Q bs)< where Q rii is the 
angular size of the black hole's influence radius and Q„i, s is the observational 
resolution. Filled symbols have r > Q i, s (r;, resolved) and open circles have 
8 r . < 9 jj (r;, unresolved). Values of T r (r/,) in the unresolved galaxies should 
be considered approximate since the luminosity profiles in these galaxies are 
not well known at r < r/, . 



M,sj3x 1O 6 M ) are the most massive to be associated with 
collisionally-relaxed nuclei, then the associated cusps could 
be resolved to a distance of ~ 0.7 Mpc, roughly the distance to 
M32 - consistent with the statement made above that the cusp 
in M32 is only barely resolved. Hence, collisional cusps are 
unlikely to be observed in galaxies beyond the Local Group. 

Unresolved density cusps might appear as pointlike nu- 
clei, particularly in dE galaxies which have low central sur- 
face brightnesses. Pointlike nuclei are in fact nearly ubiq- 
uitous in elliptical galax ies as faint as My ss —18, disap- 
pearing for My k5 — 13 (vandenBereh 1986). Luminosities 
of the nuclei are observed to average ~ 0.003 ti mes that of 
their h ost galaxies, albeit with considerable scatter (ICote et al.l 
120051) . As shown in §4, Bahcall-Wolf cusps entrain a mass 
of order 0.1M.. If the ratio of black hole mass to stellar 
mass that c haracterizes bright galaxies, M,/M ga [ w 0.0013 
( Merritt & Ferrarese 2001 ), also holds for dE galaxies, the lu- 
minosity associated with the cusps would be only ~ 10~ 4 L ga /, 
too small to explain the majority of the observed nuclei. On 
the other hand, essentially nothing is known about the masses 
(or even the existence) of black holes in spheroids fainter than 
My = — 18 (with the exception of M32, probably a special 
case) and it is possible that M. J> 10~ 3 M ?a / in these galaxies. 
It is intriguing to speculate that the disappearance of pointlike 
nuclei in dE galaxies fainter than My — 13 might signal the 
disappearance of the black holes. 

We note that the nuclear cusps of the Milky Way and M32 
extend approximately as power laws out to radii far beyond r/,. 
Even these more extended cusps would be unresolved beyond 



the Local Group and might contain enough light to explain the 
pointlike nuclei. The existence of these extended power-law 
cusps - which presumably have little to do with the presence 
of the black holes - suggests that other mechanisms, e.g. star 
formation triggered by gas infall, might be as effective as stel- 
lar dynamics at re-generating cusps. 

The black holes in the Milky Way and M32 are 
among the smallest w ith dynamically-determined masses 
(Ferrarese & Ford 2005). If smaller black holes do not exist, 
Figure [8] suggests that Bahcall-Wolf cusps might be present 
only in a small subset of spheroids containing black holes 
with masses 10 6 M Q <; M. < 3 x 10 6 M Q . However it has been 
argued that some late-type spiral s host AGN with black hole 
masses as low as ~ 10 4 M Q JHoll2004h . If so, Figure|8]sug- 
gests that Bahcall-Wolf cusps would be present around these 
black holes. 

The presence of cores, or "mass deficits," at the cen- 
ters of bright elliptical galaxie s has been taken as ev- 
idence of past merger even ts ( Milosavli evic et ail 120021 
iRavindranath. Ho & Filippenkol l20()2t lGrahamll2004l) . Mass 
deficits are o bserved to disappear in galaxies fainter than 
My w -19.5 (Milo savlievic et alJl2002t> . Could this be due 
to cusp regeneration? Figure [8] suggests an alternative expla- 
nation. Galaxies fainter than My — —19.5 are mostly unre- 
solved on scales of r/,, which is also the approximate size of 
a core created by a binary supermassive black hole. The lack 
of mass deficits in galaxies with My — 19.5 probably just 
reflects a failure to resolve the cores in these galaxies. 

5.3. Black Hole Feeding Rates 

The low-luminosity galaxies most likely to harbor Bahcall- 
Wolf cusps (Figure [8} are the same galaxies that would 
dominate the overall rate of stellar tidal disrupti ons, assum- 
ing o f course that they contain black holes (IWang & Merritt! 
2004). Published estimates of N, the rate of stellar disrup- 
tions, in such low-luminosity galaxies llSver & yim erlll999t 
Mago rrian & Tremainell999l:IWang & Merritl2004l) have al- 
most always been based on an inward extrapolation of lumi- 
nosity profiles measured at r > r/,. In principle, knowing that 
p(r) has the Bahcall-Wolf form near the black hole should 
allow a more accurate estimate of N in the low-luminosity 
galaxies that dominate the overall flaring rate. 

Here we show that the presence of a Bahcall-Wolf cusp im- 
plies a lower limit on N, of order 10~ 4 yr _1 . The stellar den- 
sity in the cusp is 



Pto 



P(n>) - 
jo 



-7/4 



(15) 



(equation |8jl, with ro = ocr/,, a w 0.2. We can write p(ro) = 
KM./rf^ where the constant K depends on the form of p(r) 
at r > ro; assuming a p ~ r~ 2 power law for r > ro, as in 
the Milky Way and many other low-luminosity spheroids, we 
find K w 4.0. The rate at which stars in the cusp are fed to the 
black hole is approximately 

4k r ar i' p , 

Ncusp w — / „,_,\> tn s **dr (16) 



m+Jo r,-ln(2/6, c .) 

(Light man & Shapiro! 1977HSver & UlmeJl999i) . Here 9, e w 
J r t /r is the angular size of the loss cone at radius r and r t 
is the tidal disruption radius, r, w (M./m*) 1 / 3 ^. This ex- 
pression assumes that the feeding rate is limited by diffusion, 
i.e. that the loss cone is "empty"; an equivalent statement 
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is that r cr it, the radius above which a star can scatter in and 
out of the loss cone in one orbital period, is greater than ocr/,. 
In the case of the Milky Way black hole, it can be shown 
that 0.2r/, < r cnf < r/,. Taking the slowly-varying logarith- 
mic terms out of the integral, we find a feeding rate for stars 
in the cusp: 

In A /GM.\ 1/2 

"-""Rwbr ) • <17) 

Evaluating 0; c at r = ocr/, and setting InA = 15, this becomes 

N cusp « 7 x lO-VXV^jS (18) 

where M,.mw and r^MW are in units of the values quoted 
above for the Milky Way. Thus, the flaring rate due to stars 
in the Milky Way cusp is ~ 10~ 4 yr . This is of course 
a lower limit on the total N since it ignores the contribution 
from stars outside the cusp, at r > uro w 0.7 pc. In fainter 
spheroids, the M. — a relation, combined with equation (II 81 . 
implies A7 « M~ 0A and hence even higher flaring rates. 

The Bahcall-Wolf solution will break down at radii where 
the physical collision time is shorter than the diffusion time 
ln(2/0; c )r, . Adopting the standard expression for the colli- 
sion time, 

T coU =[l6Vnnar^l + &)Y 1 (19) 

with © = Gm ic /(2a 2 r^) and n the number density of stars, we 
find that physical collisions begin to affect the stellar distri- 
bution at r <, 0.08 pc <, 0.023r/, for Solar-type stars in the 
Galactic nucleus. 

5.4. Gravitational Lensing 

The central parts of galaxies can act as strong gravita- 
tional lenses; the lack of a "core" image in observed lens sys- 
tems implies a lower limit on the stellar density of the lens- 
ing galaxy w ithin the central ~ 10 2 pc (iRusin & Mai 1200 fl 
lKeetonl l2003). Broken power-law density profiles like those 
in equ ation (|8j ha ve been used to mo del lensing galaxies 
(Mun oz et alJEoO 11 [Bowman et alJ2004l) . although the break 
radii in these studies were chosen to be much larger than the 
value ro w 0.2r/, that describes the Bahcall-Wolf cusps (Fig- 
ure|5}. However the presence or absence of the cusps should 
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have little effect on the lensing properties of galaxies, because 
the mass contained within the cusp is small compared with 
M t , and because even the supermassive black holes contribute 
only slightly to the lensing signal (Rusinetal. 2005). The 
low-luminosity galaxies that are likely to contain cusps (Fig- 
ure[8} are also unlikely to act as lenses. 

5.5. Dark Matter 

The distribution of dark matter on sub-parsec scales near 
the center of the Milky Way and other galaxies is relevant 
to the so-called "indirect detection" problem, in which infer- 
ences are drawn about the properties of particle dark matter 
based on measurements o f its self-annihililation by-products 
jBertone & Merritt! 12005). A recent detection of TeV ra- 
diation from the Gala ctic center by the HESS consortium 
( Aharonian et al.1 l2004l) is consistent with a particle annihila- 
tion signal, but only if the dark matter density in the inner few 
parsecs is much higher than predicted by an inward extrapola- 
tion of the standard, ACDM halo models JHoooer etaLl .2004). 
One possibility is that the dark matter forms a steep "spike" 
around the black hole (iGondolo & Silklll999l) . Particle dark 
matter would not spontaneously form a Bahcall-Wolf cusp 
since its relaxation time is extremely long. However, once a 
cusp forms in the stars, scattering of dark matter particles off 
of stars would redistribute the dark matter in phase spa ce on a 
time scale of order r r (r/,), the star-star relaxation time HVIerrittl 
2004). The ult imate result is a p ~ r~ 3 / 2 density cusp in 
the dark matter dGnedin & Primack 2004J), but with possibly 
low normalization, particularly if the dark matter distribution 
was p reviously modified by a binary black hole (tMerritt et alJ 
2002). The A/-body techniques applied here would be an ef- 
fective way to address this problem. 
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